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SUMMARY 


HEAT  SINKS:  A STUDY  OF  VARIABLE  TYPES  OF  HEAT  SINKS. 

EFFECTS  OF  VARIOUS  FACTORS  ON  TEMPERATURE  PROFILES 

U.  Fabricius 


(W 


<+ 


> 


Twelve  different  types  of  heat  sink  are  included  in  this  investigation, 
which  stresses  comparison  between  the  profiles  rather  than  comparison  be- 
tween products  of  different  manufacturers. 

The  test  plan  includes,  for  instance,  measurements  of  temperature  rise 
depending  on  load.  In  connection  with  this  measurement  the  importance  of  the 
heat  sink's  having  a black  anodized  finish  rather  than  a bright  surface  was 
examined . In  addition,  the  temperature  distribution  of  the  heat  sink  was 
studied . 


Besides  the  practical  measurements,  some  theoretical  studies  were  carried 
out,  and  some  hypotheses  out  forward  whose  correctness  was  verified  by  tests  and 
measurements.  Among  other  things,  it  was  found  that  cooling  of  the  heat  sink 
by  convection  is  strongly  dependent  on  the  space  between  the  fins. 

Computer  calculations  of  the  thermal  resistance  and  the  temperature  gra- 
dient of  the  heat  sink  were  carried  out.  Those  calculations  were  used  to  find 
out  which  proportion  between  height  of  profile  and  height  of  fin  gives  the 
lowest  thermal  resistance  compared  to  the  space  taken  up  by  the  profile. 

A simple  method  of  calculation  is  indicated  to  find  out  the  thermal  re- 
sistance of  any  cooling  profile. 

Finally,  a study  was  done  of  how  much  power  can  be  added  to  the  heat  sink 
in  intermittent  service  without  influencing  the  junction  temperature  of  the 
power  component  to  exceed  a maximum  permissible  temperature.  It  appears  that 
in  certain  cases  the  influence  tolerated  is  far  above  the  permissible  influ- 
ence according  to  the  data  sheet. 


This  study  was  carried  out  with  the  support  of  the  Technology  Council 
and  the  Society  for  Testing  and  Environmental  Engineering  (Sammenslutningen 
for  TaiiJci igheds-  og  Milj0teknik] . 
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1,  INTRODUCTION 

Thermal  problems  with  electronic  equipment  rarely  attract  attention  un- 
til too  late  in  the  development  phase.  Frequently  no  attempt  is  made  to 
solve  the  problems  until  the  pi'ototype  stage,  with  the  result  that  to  solve 
the  thermal  problems  at  all  it  is  necessary  to  adopt  a considerably  more 
costly  solution. 

We  may  hove  to  be  content,  willingly  or  unwillingly,  if  the  equipment 
works.  But  many  internal  and  external  conditions  may  cause  it  not  to  work 
satisfactorily  or  adequately;  in  stowing  the  equipment  away  it  is  necessary 
to  pay  close  attention  to  the  maximum  ambient  temperature  and  exposure  to 
solar  radiation,  and  blocking  of  radiation  bv  bpxing  the  apparatus  in. 

A poor  design  can  reduce  the  ’ifetime,  lower  the  long-term  stability, 
and  decrease  the  efficiency  of  the  equipment.  At  the  worst  it  may  result  in 
gross  defects,  actually  dangerous  products,  and  costly  or  even  fatal  damage. 

It  may  be  stated  as  a rule  of  thumb  that  the  lifetime  of  a transistor, 
for  example,  is  halved  by  a temperature  rise  of  about  1S°C. 

The  problem  of  dimensioning  and  choice  of  the  right  cooling  disk  is  a 
part  of  the  complex  system  that  is  called  thermal  design. 

What  will  the  cooling  plate  look  like  that  exactly  satisfies  the  re- 
quirements of  the  job  in  hand?  How  big  will  it  be?  Shall  it  be  bulky  or 
slender?  High  or  low?  Broad  or  narrow?  Black  or  shiny? 

If  the  model  is  in  the  prototype  stage,  the  possibility  cf  choice  is 
often  greatly  restricted,  whereas  there  is  considerably  more  freedom  if  the 
choice  is  made  at  an  earlier  stage  of  development. 

This  report  is  intended  to  cover  certain  theoretical  considerations  and 
practical  experiments  which  will  give  the  designer  a better  chance  to  choose 
the  right  cooling  plate  for  a given  purpose. 

2,  THE  TYPES  INVESTIGATED 

The  investigation  includes  ten  different  types  made  by  two  manufactur- 
ers and  two  types  of  unknown  make.  Nine  of  the  types  are  anodized  hlncV, 
one  is  shiny,  and  two  types  are  represented  by  both  a black  anodized  and 
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a shiny  type. 

In  the  selection  of  types  for  study,  the  primary  emphasis  was  placed  on 
studying  the  various  existing  designs  within  a curtain  power  range, 
rather  than  on  making  a comparative  study  of  brands. 

The  types  studied,  generally  speaking,  cover  t.he  range  within  the  field 
of  cooling  plates  which  is  used  for  cooling  low-power  components  which  give 
off  up  to  about  100  watts  as  heat,  e.g.  diodes,  transistors,  thyristors,  and 
triacs . 

These  cooling  plates  cover  the  thermal  resistance  range  from  about 
0.5°C/W  to  about  3.5°C/W. 

A list  of  variants  of  the  types  studied  is  given  in  the  table  below. 

The  thermal  resistances  shown  are  measured  values. 


2.1. 

No. 

Table  of  Types  Investigated 

Manufacturer  Type  Danish 

r Supplier 

Anodized? 

Thermal 
Resistance  at 
AT  = 100°C 

Price  per 
Hundred 
1 Sept.  1974 

1 

Fischer 

SK07 

Inotec 

Yes 

3.1 

4.93  kr. 

2 

Wakefield 

623K 

Scansupply 

Yes 

2.3 

8.90  kr. 

3 

Fischer 

SK01 

Inotec 

Yes 

2.3 

5.36  kr. 

3a 

Fischer 

SKOl 

Inotec 

No 

2.9 

4.21  kr. 

4a 

Fischer 

SK16 

Inotec 

No 

2.7 

4.79  kr. 

5 

Wakefield 

NC-303K 

Scansupply 

Yes 

2.3 

9.25  kr. 

6 

Unknown 

Yes 

1.7 

7 

Fischer 

SK02 

Inotec 

Yes 

1.3 

11.89  kr. 

7a 

Fischer 

SK02 

Inotec 

No 

1.6 

9.72  kr. 

8 

Fischer 

SKIS 

Inotec 

Yes 

1.2 

30.31  kr. 

9 

Fischer 

Kl.l 

Inotec 

Yes 

1.1 

19.72  kr. 

10 

Unknown 

Yes 

1 .1 

11 

Wakefield 

431K 

Scansuppl y 

Yes 

0.7 

17.05  kr. 

12 

Wakef ield 

465K 

Scansupp  y 

Yes 

0.6 

37.80  kr. 
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3.  TEST  PLAN 


3.1.  Survey 

The  block  diagram  below  shows  the  test  plan  used,  which  consists  in  part 
of  practical  studies  and  in  part  of  theoretical  studies. 

For  the  two  studies  in  question  --  the  dependence  of  the  convection  con- 
stant on  the  interval  between  the  ribs  and  the  intermittent  operation  of  the 
cooling  plate  --  this  means  that  certain  calculations  were  done  and  a 
hypothesis  was  set  up,  which  was  verified  by  tests. 

A more  detailed  description  of  these  studies  and  tests  can  be  found 
under  the  respective  headings. 


PRACTICAL  STUDIES 


THEORETICAL  STUDIES 


4 4.  POSSIBILITIES  OF  HEAT  TRANSFER 

lor  any  electronic  equipment,  to  which  power  is  applied,  this  power  can 
he  divided  into  Effective  power  and  waste  power-,  i.e.,  there  is  a question 
of  the  efficiency  of  the  equipment.  The  waste  power,  which  often  amounts  to 

ij  more  than  S0">  of  the  power  supplied,  must  be  removed  by  heat  transfer.  There 

li 

| 
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are  three  basic  principles  of  heat  transfer,  namely  conduction,  radiation, 
and  convection. 

A brief  description  of  their  nature  is  given  below.  If  further  in- 
formation is  desired,  the  reader  is  referred  to  ECR-15. 

4.1.  Heat  Conduction 

All  materials  allow  more  or  less  heat  transfer  by  conduction,  depending 
on  the  specific  thermal  conductivity  of  the  material. 

For  a given  material  the  amount  of  heat  transferred  by  conduction  is 
dependent  on  the  temperature  difference  as  well  as  the  dimensions  of  the 
material  and  its  specific  thermal  conductivity. 


where  R^  is  the  thermal  resistance.. 


[W] 


Starting  from  the  specific  thermal  conductivity  of  the  material,  the 
thermal  resistance  to  conduction  can  be  defined  as: 
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where  L is  the  length  of  the  material  in  m, 

A is  the  area  of  cross  section  of  the  material  in  m2  , and 

K is  the  specific  thermal  conductivity  in  W/m/',C, 

4,2.  Radiation 

Heat  transfer  by  radiation  occurs  with  the  emission  of  electromagnetic 
radiation.  For  the  temperatures  which  actual ’y  occur  in  connection  with 
electronic  equipment,  the  wavelength  falls  in  the  lower  end  of  the  infrared 
range , 


The 

heat  given  off  by  radiation  can  be  expressed  by 

the 

formula 

ps-c  • - ItSb']. 

where  C 

is  the  radiation  number 

i.n 

W/m2/k4 , 

A 

is  the  area  of  the  heat-emitting  surface 

in 

o 

m-, 

T 

is  the  temperature  of  the  object 

in 

°C,  and 

T 

is  the  radiation  temperature  of  the  environment 

in 

°c. 

6 


The  formula  is  valid  on  the  assumption  that  the  object  is  completely 
surrounded  by  the  surface  which  is  to  receive  the  heat  radiation  and  that 
that  surface  has  a high  emission  factor  and  an  area  that  is  much  greater 
than  the  area  of  the  object  in  question. 

These  conditions  will  almost  always  be  satisfied. 

4,3.  Convection 

If  a warm  body  is  placed  in  motionless  air,  the  air  is  heated,  and  the 
resultant  change  in  specific  gravity  sets  the  air  in  motion.  New  cold  air 
is  thus  brought  in,  and  a transfer  of  heat,  is  under  way.  This  mechanism  of 
heat  transfer  is  called  free  convection  or  natural  convection. 

Heat  transfer  by  free  convection  may  be  ex- 
pressed by  the  following  formula: 

PK  = K • C • A • HZ.  • AT. 

where  K is  a constant, 

C is  a constant  determined  by  the  orienta- 
tion and  geometric  shape  of  the  object, 

A is  the  area  of  the  heat-emitting  surface 
in  m-0  , 

AT  is  the  difference  in  temperature  between 
the  heat -emitting  surface  arid  the  object 
in  °C,  and 

L is  the  characteristic  dimension  in  meters. 

f,.  TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

On  the  pages  that  follow,  the  rise  in  tempera- 
ture is  shown  as  a function  of  the  power  applied. 
All  measurements  are  done  at  an  ambient  temperature 
of  20°C  under  free-field  conditions;  i.e.,  in  the  measurement  chamber  the 
only  circulation  of  air  is  that  which  is  created  by  free  convection  from  the 
cooling  plate. 

During  the  experiment  the  cooling  plate  was  freely  suspended  in  the 
room,  with  such  a distance  between  the  various  subjects  of  the  experiment 


Figure  1.  Isotherms 
Around  a Horizontal 
Pipe  With  Free  Con- 
vect ion . 


that  they  could  not,  from  the  thermal  point  of  view,  "see"  each  other. 

The  individual  measurement  points  are  shown  on  the  sketches  which  ac- 
company the  graphs.  The  temperature  was  measured  with  iron-constantan 
thermocc'Uples . 

For  the  various  types  curves  are  presented  for  the  relatively  warmest 
point  and  coldest  point  measured  on  the  cooling  plates.  Where  the  type  is 
under  study  in  both  a bright  and  a black  anodized  version,  the  two  sets 
of  curves  are  drawn  on  the  same  graph  to  facilitate  comparison.  In  addi- 
tion, the  curve  given  in  the  data  sheet  is  drawn  in  if  available. 

5.1.  Comparison  With  the  Specifications  of  the  Data  Sheet 

For  most  types  it  is  found  that  the  measured  values  fall  within  those 
given  in  the  data  sheet;  e.g.,  No.  1 and  No.  2 are  considerably  better 
than  their  respective  data  sheets  indicate.  It  should  be  noted  that  these 
two  are  small  cooling  plates.  No.  8 is  poorer  than  its  specifications. 

No.  9,  which  in  the  data  sheet  is  merely  specified  with  AT  = 100°C,  meets 
the  specification  at  that  point,  but  below  AT  = loO°C,  where  the  cooling 
plate  would  normally  be  used,  the  thermal  resistance  is  higher  than  stated 
(Ky  = 1.1°C/W),  No.  8 and  No.  9 are  relatively  large  cooling  plates. 

There  is  thus  a tendency  for  small  cooling  plates  to  be  better  than  is 
claimed  for  them,  while  larger  cooling  plates  are  poorer  than  claimed. 

5.2.  Comparison  Between  Bright  and  Black  Anodized  Finish 

Comparison  of  the  bright  and  the  anodized  versions  shows  a rise  in 
thermal  resistance  of  about  20  to  2500  for  the  bright  versions  as  against 
the  anodized  ones. 

Since  the  surface  treatment  can  have  significance  only  for  heat 
transfer  by  radiation  --  the  emission  factor  is  reduced  from  about  5.5 
(anodized)  to  0.5  (bright)  --  we  can  get  an  impression  from  this  of  the 
part  that  radiation  plays  in  the  heat  transfer  for  the  types  available  in 
both  an  anodized  and  a bright  version. 

5,5.  Temperature  Cradicnts  in  the  Heat  Sink 

The  table  below  shows  the  maximum  measured  temperature  gradient  in 
the  heat  sink  itself  at  a maximum  absolute  temperature  of  100°C,  The  less 


H 


the  temperature  gradient,  the  more  fully  the  cooling  plate  is  utilized  and 
the  more  effective  it  is. 


Type  No.  1 2 3 3A  4A  5 6 7 7A  8 9 10  11  12 

Temperature  Gradient  in  °C  9 14  14  11  8 21  13  / 7 37  26  31*  18  20* 

*At  80°C , 


The  next  section,  on  thermovision,  gives  further  information  about 
the  temperature  distribution  on  the  various  cooling  plates. 


6.  THERMOVISION 


The  pages  following  the  temperature  gr-  ent  graphs  show  thermovision 
pictures  of  the  various  types.  It  should  be  noted  that  No.  4 and  No.  6 
are  not  included. 


The  thermovision  pictures  were 
taken  with  the  AGA  thermovision  sys- 
tem 680.  That  apparatus  operates  at 
a wavelength  of  2 to  5.6  ym,  and 
without  filters  it  can  measure  tempe- 
ratures from  -30°  to  +850°C.  The 
apparatus  scans  with  a frequency  of 
16  pictures  per  second. 

In  the  present  case  we  are  not 
concerned  with  measuring  the  absolute 
temperature,  but  merely  with  getting 
an  impression  of  the  temperature  dis- 
tribution on  the  cooling  plate. 


Figure  2.  AGA  Thermovision 

Apparatus,  Type  680.  It  should  be  mentioned  that  the 

scale  is  not  the  same  in  the  various 
pictures;  the  object  was  simply  focused  on  in  such  a way  as  to  fill  out 
the  frame  as  well  as  possible. 


Since  the  apparatus  measures  the  radiation  of  heat,  it  is  very  diffi- 
cult to  make  measurements  on  the  bright  cooling  plates,  which  have  a very 
low  emission  factor  (0.2  to  0.5)  and  a very  high  reflection  factor  (about 
0.95).  The  measured  value  is  the  sum  of  the  characteristic  radiation  and 
the  reflection;  see  Figure  3 [page  22] . 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

Frse-FIeld  Measurements 


Number  ! 

Manufacturer:  Fischer 

Type:  SK  07 


Ambient  Temperature  ■ 20 
Air  Velocity  = 0 m/second 
Volume  « 76.1  cm3 


20 


TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

Free-FIeld  Measurements 

No.  2 Ambient  Temperature  - 20°C 

Manufacturer:  Wakefield  Air  Velocity  =■  0 m/second 

Type:  623  Volume  * 98. d6  cm3 


M 


$ 


1 


220 

200 

180 


\c/\ 


140 


120 


100 


80 


eo 


Rf t*s  Parallel 
to  Vertical 


ft 


"■a  r&i 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

Free-Field  Measurements 

No.  6 Ambient  Temperature  = 20°C 

Manufacturer:  unknown  Air  Velocity  ■ 0 m/second 

Type:  unknown  Volume  - 297.6  cm3 


5j 

r: 

t 

it: 


f 

% 


i 


H 


\ 

k BLACK  ANODIZED 
I SURFACE 

J 


Ribs  Paral lei 
to  Vertical 


■V  r 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 


No.  7 

Manufacturer:  Fischer 

Type:  SK  02 


Free-Fleld  Measurements 

Ambient  Temperature  ■ 20°C 
Air  Velocity  «•  0 m/second 
Volume  - 43*1.7  cm 3 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 


Free-Field 

No.  9 

Manufacturer:  Fischer 

Type:  K I . 1 


Measurements 

Ambient  Temperature  * 20°C 
Air  Velocity  * 0 m/second 
Specification:  1.1°C/W 

Volume  * 6D?.$  cm3 


[*c] 


o 


/ 


© 

© 


> 


BLACK  ANODIZED 
SURFACE 


DATA  SHEET  SPECIFICATION 
BLACK  ANODIZED  SURFACE 


Ribs  Para  1 le) 
to  Vertical 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 


Free-Field  Measurements 


No.  10 

Manufacturer:  unknown 

Type:  unknown 


Ambient  Temperature  * 20 
Air  Velocity  = 0 m/second 
Volume  ■ 633.7  cm3 


-T  [«C] 
200  ■ 


Measuring  Points 

© 9 ® 


mXl  t.  20 


rr 


Height  6.9  cm 
Measured  With 
Ri bs  Paral lei 
to  Vertical 


BLACK  ANODIZED  SURFACE 


20  40  60  OO  100  120  140  P [tf] 


TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

Free-Field  Measurements 


No.  1 1 

Manufacturer: 
Type:  A3 1 A 


Wakefield 


Ambient  Temperature  « 20°C 
Air  Velocity  **  0 m/second 
Volume  - 695  cm3 


BLACK  ANODIZED 
SURFACE 
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TEMPERATURE  RISE  AS  A FUNCTION  OF  LOAD 

Free-Field  Measurements 

Ambient  Temperature  = 20°C 
Air  Velocity  ■ 0 m/second 
Volume  * 1,311  cm3 

0 


No.  12 

Manufacturer:  Wakefield 

Type:  A65 


T,>T2 


MEASURED 


Tigure  3.  In  Measurements  of  Thermal 
Radiation,  What  is  Measured  is  Charac- 
teristic Radiation  + Reflection, 


It  is  evident  that  in  the 
main  the  small  types  assume 
the  same  temperature  throughout 
their  surface,  while  the  tem- 
perature gradient  on  the  larger 
ones  is  more  pronounced,  espe- 
cially in  the  case  of  No.  8 and 
No,  lo.  This  is  in  harmony 
with  the  difference  measured 
under  5.4  above. 


Fischer  type  SK07 


Wakefield  type  623 


1 


Fischer  type  SKI'  1 


Waketield  type  s(H 


7.  COMPARISON  OF  QUALITY 

A comparison  of  quality  may  be  based  on  the  principle  thpt  the  least 
possible  thermal  resistance  is  desired  within  a given  space. 

If  we  plot  the  thermal  resistance  of  the  various  types  as  a function 
of  their  volume  on  a simple  logarithmic  coordinate  system,  it  turns  out 
that  we  get  a s*  aight  line.  What  is  meant  here  by  tl^e  volume  of  the  heat 
sink  is  the  volume  or  the  space  that  the  heat  sink  takes  up. 


The  table  below  gives  the  characteristic  values  for  the  various  types. 


No . 

Manu- 

facturer 

Thermal  Re- 
Type  sistance  at 
AT  - 100°C 

Thermal 

Capac- 

ity 

Thermal 

Time  Volume  Weight 
Const , 

Ano- 

dized? 

"Specific 

Gravity" 

cC/W 

Ws/°C 

min , 

cm3 

grams 

g/cm? 

1 

Fischer 

SK07 

3.1 

80 

4.2 

76 

89 

yes 

1.17 

2 

Wakefield 

623 

2.3 

98.5 

3.8 

99 

110 

yes 

1.11 

3 

Fischer 

SK01 

2.3 

111 

4.2 

181 

124 

yes 

0.68 

3a 

Fische. 

SK01 

2.9 

111 

5.4 

181 

124 

no 

0.68 

4a 

F ischer 

SK16 

2.7 

140 

6.4 

200 

156 

no 

0.78 

5 

Wakefield 

303 

2.3 

146 

5.6 

198 

163 

yes 

0.82 

6 

Unknown 

1.7 

183 

5.1 

2S8 

205 

yes 

0.69 

7 

Fischer 

SK02 

1 .3 

244 

5.4 

435 

272 

yes 

0.63 

7a 

Fischer 

SK02 

1.6 

244 

6.5 

4 35 

272 

no 

0.63 

8 

Fischer 

SKI  5 

1.2 

753 

14.6 

675 

841 

yes 

1.25 

9 

Fischer 

K1 .1 

1.1 

S65 

10.4 

608 

631 

yes 

1 .04 

10 

Unknown 

1.1 

269 

5.0 

634 

300 

yes 

0.47 

11 

Wakefield 

431A 

0.6 

334 

4 .8 

695 

373 

yes 

0.54 

12 

Wakefield 

465 

0,6 

722 

7.6 

1 ,311 

806 

yes 

0.61 

The  straight  line  (see  Figure  4)  can  be  used  as  a basis  for  t compara- 
tive evaluation  of  quality,  since  we  can  say  that  for  heat  sinks  whose 
measurements  fell  above  the  straight  lino  the  thermal  resistance  is  great, 
while  for  those  which  fail  below  the  line  the  thermal  resistance  is  low, 
everything  being  taken  in  relaticr  to  the  space  taken  up  by  the  heat  sink. 


?.k 


m 


battik 


attmiAimm 


At  100*c 


O Specified  Value 
+ Measured  Value 
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Figure  4.  Thermal  Resistance  as  a Function  of  the  Volume. 

For  comparison,  the  thermal  resistance  is  plotted  in  for  a black  ano- 
dized aluminum  cube  and  also  for  a similarly  black  anodized  A1  plate,  wher< 
the  thermal  resistance  from  the  center  of  the  plate  to  the  outer  edge  is 
only  one  third  of  the  thermal  resistance  from  the  plate  to  the  surrounding 
air.  These  two  curves  may  be  considered  as  extreme  limits,  since  the  cube 
is  one  of  the  worst  possible  designs,  from  the  thermal  point  of  view,  and 


the  plate  is  the  design  that  may  be  considered  in  practice  as  the  best. 

There  was  no  attempt  to  optimize  the  thickness  chosen,  however. 

For  small  values  of  thermal  resistance  the  area  of  the  plate  is  im- 
practically  large,  so  that  it  does  not  represent  a realistic  solution.  It 
is  included  here  only  in  order  to  provide  an  extreme  limit. 

A general  appraisal  shows  that  small  heat  sinks  lie  in  the  range  be- 
tween a plate  and  a cube,  approaching  the  cube  for  a size  of  about  1°C/W  in 
thermal  resistance  and  then  shifting  back  toward  a plate  design.  If  we  look 
at  the  geometric  form  of  the  heat  sinks,  too,  we  see  that  small  heat  sinks 
approach  a plate  in  design,  while  the  medium  size  with  a thermal  resistance 
of  about  1°C/W  is  more  massive  and  approaches  a cube  in  shape.  Lastly,  the 
large  sizes  (No.  11  and  No.  12)  have  a design  in  which  the  individual  cool- 
ing fins  are  set  at  such  a distance  from  each  other  and  in  such  a geometric 
arrangement  that  they  may  be  regarded  as  composite  plates. 

In  conclusion  it  should  be  noted  that  the  two  limiting  curves  are 
purely  theoretically  calculated  curves,  based  on  the  assumption  that  the 
whole  mass  has  the  same  temperature. 

8.  THEORETICAL  INVESTIGATIONS  OF  THE  HEAT  SINK 

The  preceding  sections  are  based  on  practic  :1  measurements.  In  the 
following  sections  theoretical  reflections  and  investigations  will  be  under- 
taken, both  in  order  to  facilitate  the  calculation  of  a heat  sink  and  in 
order  t"  approach  a possible  optimal  design  of  a heat  sink. 

8.1.  The  Convection  Resistance  as  a Function  of  the  Distance 

The  heat  released  from  a warm  body  by  convection  may  be  expressed  by 
the  formula 

P = K • C • A • 4/ljl  • AT. 

This  expression  contains  a number  of  parameters,  some  of  which  are  mea- 
surable and  others  very  difficult  to  measure. 

It  is  generally  true  for  electronic  apparatus,  however,  that  the  "con- 
stant" K,  which  covers  many  of  the  quantities  that  are  hard  to  measure,  can 
be  set  equal  to  2.5,  so  that  we  get  the  simplified  expression 
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4/37 

P = 2.5*C*a/^-*AT. 

However,  for  heat  sinks,  among  other  things,  we  cannot  always  allow 
ourselves  to  simplify  the  expression  as  shown  above.  It  is  found  that  when 
we  have  two  heat-emitting  surfaces  parallel  to  each  other,  the  factor  K is 
gre.-*ly  dependent  on  the  distance  between  the  two  surfaces. 

If  we  look  at  the  air  velocity  curve  for  a vertical  heated  surface  with 
free  convection  (Figure  5) , we  see  that  the  curve  reaches  a maximum  value 
at  a distance  of  about  l/rf  from  the  heat-emitting  surface,  and  the  velocity 
drops  to  zero  at  the  distance  y = 3. 


Figure  5.  Temperature  and  Air  Velocity  Plotted  as 
Functions  of  Distance  From  a Heat -Emitting  Surface 
With  Free  Convection. 


This  limiting  distance  9 for  a laminar  flow  is  defined  by  the  follow- 
ing expression: 

3lam  = 5.83  / 


/rTT~T 


where  9lam 

i s 

the 

limiting  distance  in 

laminar  flow 

1 

m 

6 

is 

the 

dynamic  viscosity  of 

the  air 

1 

s 

1 

is 

the 

distance  to  the  edge 

of  inflow 

1 

m 

v 

is 

the 

velocity  of  the  airflow 

1 

m 

s* 
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Figure  6.  Graphic  Representation  of 
the  Limiting  Distance  as  a Function 
of  the  Height  and  Temperature  of  the 
Surface . 


Now  if  we  consider  the  air 
velocity  between  two  heat-emit- 
ting surfaces,  we  will  theoret- 
ically get  curves  like  those 
shown  in  Figure  7,  We  shall  now 
show  that  in  practice  we  cannot 
allow  an  overlapping  of  the  two 
velocity  curves,  if  free  convec- 
tion is  to  be  maintained. 

To  establish  the  exact  de- 
pendent relationship  between  the 
heat  that  can  be  given  off  in 
free  convection,  i.e.  the  convec- 
tion resistance,  and  the  distance 
between  two  heat -emitting  sur- 
faces, experiments  were  carried 
out  in  which  a constant  and  equal 
amount  of  heat  was  applied  to  two 
parallel  surfaces  and  the  tempera- 
ture rise  was  measured. 

Figure  8 is  a sketch  of  the 
experimental  arrangement.  By  ap- 
plying equal  amounts  of  heat  to 
the  two  plates,  so  that  they 
took  on  the  same  temperature,  we 
eliminated  the  possibility  of  a 
heat  transport  between  the  two 
surfaces  by  radiation.  Loss  of 
heat  from  the  back  of  the  plates 
was  reduced  as  much  as  possible 

by  effective  insulation;  emission  of  heat  from  the  back  cannot  be  entirely 
eliminated,  but  since  the  thermal  resistance  here  is  constant,  it  can  be 
eliminated  in  the  calculations. 


Figure  7,  Curve  of  Velocity  of  Move- 
ment of  Air  Between  Two  Vertical 
Heated  Surfaces. 
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POWER  SUPPLIED  TO  EACH  PLATE  P = 4 W 
A = (2.5,  1.5,  1.0,  0.8,  0.7,  0.6,  0.5,  0.4) 
Area  of  the  Plate:  A =*  30  cm2 

CONVECTION  RESISTANCE  (K  - 2.5) 

Rc  = 43. 1°C/W 

Figure  8.  Sketch  of  the  Experimental  Set-Up 
for  Study  of  the  Convection  Constant's  De- 
pendence on  the  Distance  Between  Heat- 
Emitting  Surfaces. 


8.2.  Thermal/Electric  Analogy 


The  result  of  the  ex- 
periment can  be  seen  from 
the  two  curves  of  Figure  9, 
showing  the  temperature 
rise  on  the  plates  and  the 
size  of  the  constant  K as 
a function  of  the  distance 
between  the  surfaces. 

As  the  figure  shows, 
the  temperature  begins  to 
rise  when  the  distance  be- 
tween the  plates  is  less 
than  2x9;  this  exactly 
corresponds  to  the  point 
where  an  overlapping  of 
the  air  velocity  curves  of 
the  two  surfaces  begins. 

When  the  rise  in  tempera- 
ture is  inserted  in  the 
formula  for  free  convec- 
tion, we  get  the  indicated 
dependence  of  the  size  of 
K on  the  distance  between 
the  surfaces. 


Since  part  of  the  tools  for  computing  thermal  designs  are  based  or. 
electrical  magnitudes,  and  since  these  magnitudes  are  generally  easier  to 
compute  with,  it  will  be  a great  advantage  if  we  can  convert  the  thermal 
magnitudes  into  analogous  electrical  magnitudes. 


It  turns  out  that  a very  simple  analogy  can  be  found  between  the 
thermal  and  electrical  magnitudes. 
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Figure  9.  The  Curves  Show  the  Dependence  of  the  Temperature  [T]  and  Convec- 
tion [C]  on  the  Distance  Between  Two  Heat-Emitting  Surfaces. 

The  correspondence  is  as  follows: 

Thermal  Magnitude  Electrical  Analogy 


p 

power 

in 

watts 

'\y 

I 

current 

in 

amperes 

T 

temperature 

in 

°C 

V 

voltage 

in 

volts 

R 

thermal  resistance 

in 

°C/W 

R 

el ectrical 

resistance 

in 

ohms 

C 

thermal  capacity 

in 

Ws/°C 

C 

electrical 

capacity 

in 

farads 

Transfer  of  heat  by  conduction: 


K*A 


*L  “ K • A 
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Transfer  of  heat  by  convection: 


PK  = K*  C*  A 


Transfer  of  heat  by  radiation: 


*K 


1_ 

K*C*  A 


C*  A* 


^ (Too)  ~ (Too>  ^ 


Rs  = 


,T  T * 

< > ~ iL 


* A* 1 'ToV  ~ 1 


100' 


8.3.  A Sijnpie  Calculation  of  the  Thermal  Resistance  of  the  Heat  Sink 


A very  simple  method  is  given  be- 
low for  calculating  the  thermal  re- 
sistance of  a heat  sink. 


It  is  assumed  that  the  heat  sink 
has  the  same  temperature  throughout. 
No.  7 is  used  as  an  example  for  the 
calculation . 


Figure  10.  Radiation  Area  of  No.  9, 


For  this  type,  the  area  effec- 
tive for  radiation  is  limited  to  the 
area  shaded  in  Figure  10.  That  is, 
the  area  is 


A = 2 x 0.115  x 0.075  + 2 x 0.063  x 0.075  = 0.0267  m2 . 

On  the  basis  of  Section  8.1,  we  find  the  radiation  resistance  to  be 

1 


= 


)/(T 


t2) 


Since  we  are  calculating  with  a black  anodized  surface,  C - 5.7,  and 
if  we  wish  to  find  the  thermal  resistance  at  AT  = 100  and  we  assume  an 
ambient  temperature  T;,  = 20°C?>we  get  Tj  = 120°C.  With  those  values  substi 
tuted  in  the  formula,  we  get 


100 


Rs  = 5.7  x~ 0.0267  164.9 


= 3.98°C/W. 
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In  convection  it  is  the  entire  area  of  the  heat  sink  that  contributes. 

A relatively  simple  way  of  calculating  the  area  is  as  follows  (see 
Figure  11): 

A = 2 x A2  + 16  + Aj  - 30  x A3 

= 2 x 0.075  x 0.092  + 16  * 0.075  * 0.063  - 30  x 0.075  * 0.003  = 0.0827  m2 . 


A 3 


Figure  11.  Areal  Contribution  to 
the  Area  Effective  for  Convection 
for  No.  7. 


According  to  Section  8.1,  the 
thermal  convection  resistance  is 


■'iv > 

i+Z’F 

K • C • A • 7 f 

where  the  characteristic  length  in 
this  case  is  L = 0.075  ai  and  the  con- 
stant C = 0.55,  since  the  surfaces 
are  vertical.  From  Section  8,2, 

K = 1.95;  therefore 


Rjc  - 


1 


1.95  • 0.55  • 0.0827 


lV/Too 


= 1 , 87°C/W , 


Radiation  resistance  and  convection  resistance  may  be  regarded  as  a 
parallel  relationship,  and  we  get 


KK  ’ Ks  _ 1,87  • 3,98 
rK~+  Rs  it87  + 3.98 


1 . 3°C/W. 


The  measured  thermal  resistance  at  AT  = 100°C  is  R = 1.33°C/W.  Thus 
in  this  case  there  is  close  agreement  between  the  calculated  and  the  meas- 
ured results. 

This  method  of  measurement  does  not  show  the  temperature  gradient  in 
the  heat  sink, 

The  next  section  will  show  the  application  of  a computer  program  to 
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calculation  of  the  thermal  resistance  of  the  heat  sink.  This  calculation 
gives  not  only  the  thermal  resistance  (indirectly)  but  also  the  temperature 
distribution  in  the  heat  sink. 

8.4.  Computer  Calculation 

The  principle  of  the  use  of  a computer  in  thermal  design  of  electronic 
apparatus  will  be  described  in  greater  detail  in  this  section  and  will  be 
illustrated  both  by  a little  example  and  by  calculation  of  the  temperature 
distribution  in  a heat  sink. 

The  technique  is  based  on  the  use  of  circuit  programs;  i.e.,  the  con- 
verted magnitudes  described  in  Section  8.1  will  be  used. 

Since  the  calculation  is  rather  complex  and  may  not  seem  very  clear 
at  first  glance,  the  principle  will  first  be  explained  by  means  of  a little 
hypothetical  example. 

8.4.1.  The  Principle  of  the  Use  of  Circuitry  Programs 

figure  12a  shows  a simple  example  of  a thermal  model  and  Figure  12b 
the  analogous  electrical  model. 

The  power  is  generated,  for  example,  by  a power  transistor  in  a TO-3 
housing.  The  transistor  is  mounted  in  one  end  of  the  structure  and  the 
heat  is  carried  by  conduction  down  through  the  structure,  part  of  it  being 
given  off  to  the  surroundings  by  radiation  and  convection  on  the  way. 

The  analogous  electrical  model  in  Figure  12b  is  divided  into  nodes  in 
such  a way  that  between  two  nodes  there  is  either  a resistance  or  a current 
generator.  Node  0 is  the  "frame." 

Between  nodes  1 and  2 the  resistance  corresponds  to  the  thermal  re- 
sistance from  the  junction  to  the  housing  of  the  transistor;  the  resistance 
between  nodes  2 and  3 corresponds  to  the  thermal  contact  resistance  from 
the  housing;  and  the  resistance  between  3 and  4 corresponds  to  the  thermal 
conduction  resistance  in  the  first  half  of  Section  1.  The  two  resistances 
between  4 and  0 correspond  to  the  convection  resistance  and  radiation  re- 
sistance to  the  surroundings  respectively,  and  so  on. 
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Figure  12.  Example  of  a thermal  model  (a)  and  the  analogous  electrical  model 
corresponding  to  it  (b)  . 

The  voltage  in  relation  to  the  frame,  i.e,  the  temperature  in  relation 
to  the  ambient  temperature  for  the  various  nodes  can  be  calculated  by  means 
of  a circuitry  program. 

It  will  be  found  that  the  more  sections  the  structure  is  divided  into, 
the  more  accurate  the  calculations  will  be, 

8.4.2.  Computer  Calculation  of  a Heat  Sink. 

A corresponding  description  of  the  thermal  calculation  of  a heat  sink 
will  be  given  below. 


3^ 
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Figure  13.  The  Type  of  Heat  Sink  Used 
for  the  Computer  Calculations. 

resistance  and  the  convection  resistanc 

Radiation  Resistance.  --  There  is 
tween  the  ribs,  and  since  the  radiation 
this  will  be  taken  into  account. 


Because  of  the  symmetrical 
design  of  the  heat  sink  (see  Fig- 
ure 13),  we  can  content  ourselves 
with  considering  Vg  of  the  whole, 
although  this  means  leaving  the 
difference  in  temperature  between 
top  and  bottom  out  of  considera- 
tion. 

The  heat  sink  is  divided  up 
into  sections,  such  that  the  first 
section  includes  the  part  from  the 
center  line  to  the  first  rib  in- 
clusive, the  second  section  in- 
cludes the  part  between  the  first 
and  second  ribs  plus  the  second 
rib,  etc.  The  individual  sectioi  s 
are  then  further  divided  into  a 
number  of  nodes,  and  the  necessary 
resistance  is  plotted. 

In  calculating  the  individual 
resistances,  two  "types"  will  be 
considered,  namely  the  radiation 
between  the  ribs. 

very  low  temperature  gradient  be- 
resistance  is  temperature-dependent, 


Convection  Resistance.  --  It  follows  from  what  was  said  in  Section  8.J 
that  there  is  a very  strong  dependent:  relationship  between  the  interval  be- 
tween the  ribs  and  the  convection  constant.  To  get  as  flexible  a method  of 
calculation  as  possible,  a program  was  devised  that  makes  these  calculations 

The  program  can  generate  the  necessary  input  data  for  the  circuitry 
program,  but  the  input  program  can  only  cope  with  a heat  sink  design  such  as 
is  shown  in  Figure  13. 


9k, 


' mv  -w  * wiwtwpppji^  !■  jlp  u,wpi 


Figure  14.  Analogous  Electrical  Model  of  a Heat  Sink.  Rp  = Conduction  Re- 
sistance, Rk  = Convection  Resistance,  and  Rg  = Radiation  Resistance. 


The  necessary  input  data  will  be  given  as  answers  to  the  questions 
asked  by  the  program,  such  for  example  as  the  dimensions,  radiation  number, 
and  load. 

As  a test  of  the  computer  calculation  just  described,  a calculation 
was  done  on  heat  sink  No.  7,  the  result  of  which  can  be  seen  in  Figure  16. 
As  the  graph  shows,  there  is  close  agreement  between  the  measured  and  the 
calculated  values.  At  the  coldest  measuring  point  (Point  f>)  there  is  a 
slight  discrepancy,  which  is  due  to  the  tact  that  the  outermost  rib  ot  the 
heat  sink  has  two  projections  provided  for  mounting  the  heat  sink.  Ihese 
two  projections,  which  were  not  taken  into  account  in  the  calculation, 
will  cool  the  outermost  rib  still  further. 

A more  detailed  description  of  the  input  program  is  to  be  found  in 
Appendix  1 , 
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WHAT  15  THE  PASSWORD  FOR  FILE  DCFIL. 
PASSWORD 


DETTE  ER  ET  PROCRAM  TIL  BE REG XING  AF 
TEMP=FORDELINGEN  PAA  ET  KOLEPROFIL 

HA’ AD  ER  PROFILENS  LAENGDE  - L. RIBBEHOJDEN  ~ H. 
RIBBEAFSTAXDEN  - A.AFSTANDEN  FRA  CENTERLINIEX 
TIL  FORSTE  RIBBE  - C .TYKKELSEN  AF  RIBBEN  1.  VED 
FODFA*  - 51,  OG  2.  YDERST  - S2,  MAAL  I METER.  : 

? 0.037 5* 0.03^0.007. 0.01 8. 0.003. O.ool 5 

HVOR  MANGE  RIBBER  ER  DER  PAA  HVER  HALVDEL  AF  PROFILEN 
?.  4 

HVAD  ER  KROPPENS  HALVE  TYKKFL5E 

MLLLEM  CENTER  OG  1.  RIBBE,  MELLEM  1.  OG  2.  O.S.V. 
VERDIERXE  INDTASTES  EX  AF  GANGEN 
? 0.002 
? 0.002 
? 0.002 
? 0.0012 

LAXGDE  ER  DELT  OP  I 1 KNUDEPUNKTER 

RIBBEN  ER  DELT  OP  I 3 KNUDEPUNKTER 

ER  DETTE  OK  JA  ELLER  NEJ?  NEJ 

HVAD  SKAL  ANTALLET  AF  KNUDEPUNKTER  VARE  I 

HEXHOLPSVTS  1.AENGDE  OG  RIBBF?  4,3 

HVAD  ER  OMGIVELSESTEMPERATUREN . 

HVAD  ER  KOLEPROFILENS  OVERTEMPERATUR , 

SAMT  HVAD  ER  TEMPERATURFORSKELLEN  MELLEM  RIBBERNE. 

7 29  35  1 

ER  KOLEPROFILEN  MONTERET  LODRET  ELLER  VANDRET?  LODRET 
HVAD  ER  STRAALINGSTALLET?  5.7 

HVAD  ER  VARMELEDNINGSEVNEN  FOR  MATERIALET?  200 
HVOR  STOP  ER  EFFEK^EN  I WATT  (l/8)?  3 24 


Figure  15.  The  Questions  Asked  by  the  input  program. 


8.5.  Optimization  of  the  Heat  Sink 

In  the  choice  of  a heat  sink  there  are  a number  of  factors  connected 
with  the  apparatus  in  question  which  restrict  the  possibilities  of  choice. 
The  optimal  solution  is  an  infinitely  thin  and  infinitely  conductive  plate 
with  a high  emission  factor.  This  theoretical  solution  is  obviously  impos- 
sible in  practice.  One  of  the  solutions  that  come  closest,  as  mentioned 
above,  is  a black  anodized  aluminum  plate  with  a reasonably  low  conduction 
resistance.  This  solution  is  applicable  for  the  great  majority  of  appara- 
tus, however,  only  if  a very  small  cooling  surface  is  to  be  used,  since 
otherwise  the  dimensions  of  the  plate  will  be  of  an  order  of  magnitude  that 
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Figure  16.  The  Result  of  the  Computer  Calculations. 


is  not  realistic. 

If  in  the  various  heat  sink  designs  the  plate  is  cut  up  into  larger 
and  smaller  pieces  and  put  together  in  a different  way,  an  attempt  is  made 
to  preserve  the  thermal  properties  of  the  plate  by  providing  possibilities 
of  heat  transfer  by  radiation  and  convection  while  getting  a handier  geo- 
metric pattern. 

In  this  change  from  plate  tc  profile  the  aim  must  be  the  best  solution; 
i.e..  the  greatest  possible  part  of  the  total  area  of  the  cooling  profile 
must  be  used  as  a radiation  area  (this  obviously  can  be  significant  only 
if  the  profile  is  black  anodized)  and  the  distance  between  the  ribs  must  be 
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Figure  17.  Significance  of  the  Height  of  the  Profile  and  the  Height  of  the 
Ribs . 


great  enough  so  that  the  convection  resistance  is  not  significantly  in- 
creased, but  not  so  great  as  to  result  in  poor  utilization  of  the  volume  of 
the  heat  sink.  In  addition,  the  "plates"  must  be  of  such  thickness  that 
Che  temperature  gradient  through  the  profile  does  not  become  too  great. 

What  design  and  what  mutual  ratio  between  the  dimensions  of  the  heat 
sink  will  give  the  lowest  thermal  resistance  in  relation  to  the  volume  of 
the  heat  sink,  or  in  other  words,  the  best  utilization  of  the  given  space? 

Figure  17  shows  some  calculated  curves  of  a heat  sink,  which  is  shown 
above  the  curves.  The  calculations  were  done  as  indicated  in  Section  8.4. 
The  ratio  of  the  height  of  the  profile  h to  the  rib  height  r is  taken  as  a 
parameter.  The  greater  the  ratio  the  lower  the  thermal  resistance.  This 
is  most  marked  for  small  heat  sinks  (small  volume),  while  for  the  larger 
heat  sinks  the  ratio  of  the  height  of  the  profile  to  the  height  of  the  rib 
is  of  less  significance. 

If  a star-shaped  design  is  used  instead,  as  shown  for  example  in  Fig- 
ure 18,  it  is  found  that  a ratio  of  2.5  between  the  neight  of  the  profile 
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Figure  18,  Example 
of  a Star-Shaped 
Design, 


and  the  radius  gives  a curve  that  coincides  with 
the  curve  in  Figure  17  with  the  ratio  5, 

It  may  thus  be  concluded  that,  especially  in 
the  case  of  small  profiles,  as  high  a profile  as 
possible  may  be  used  and  by  way  of  compensation  a 
smaller  rib,  while  for  larger  profiles  the  ratio 
between  the  height  of  the  profile  and  the  height 
of  the  rib  is  of  less  significance. 

It  may  also  be  concluded  that  the  star-shaped 
profile  is  a better  design  than  the  profile  with 
perpendicular  ribs. 


9.  INTERMITTENT  SERVICE 


Up  to  this  point  the  heat  sink  has  been  studied  only  where  there  was  a 
constant  input  of  heat  from  the  semiconductor,  but  often  the  latter  is  used 
in  intermittent  service. 


TJ  * .1  -hs  AT  hs 


Figure  19.  Equivalent  Diagram  for  the  Heat  Sink  and 
the  Load  to  Which  the  Heat  Sink  is  Subjected. 

This  means  that  the  heat  sink  has  a chance  to  cool  off  more  or  less 
between  the  individual  pulses,  so  that  the  same  maximum  temperature  is  not 
attained  in  the  heat  sink  and  the  semiconductor  as  under  a constant  load. 

It  is  thus  possible  to  apply  a greater  amount  of  power  in  using  the 
semiconductor  in  pulse  operation  than  in  constant  operation,  without  ex- 
ceeding a specified  maximum  temperature  at  the  junction. 
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We  shall  discuss  below  how  much  greater  load  it  is  permissible  to  put 
on  the  heat  sink  in  pulse  service  as  compared  to  steady  service. 

From  the  thermal/ electrical  equivalence,  the  analogous  electronic 
equivalent  diagram  for  the  heat  sink  can  be  drawn  [Figure  19] , where 


P is  the  power  developed  in  the  semiconductor  W 

ATj  the  temperature  rise  at  the  junction  in  com- 
parison to  the  ambient  temperature  °C 

AThs  the  maximum  temperature  rise  in  the  heat  sink 

in  comparison  to  the  ambient  temperature  °C 

Tamb  the  ambient  temperature  °C 

Rj_hs  the  thermal  resistance  from  junction  to  heat  sink  °C/W 

Rhs-amb  the  thermal  resistance  from  heat  sink  to  ambience  °C/W 

C'hs  the  thermal  capacity  of  the  heat  sink  J/°C 

t the  time  min, 

tj  the  pulse  width  or  "on"  time  min. 

T the  period  time  min. 


In  the  equivalent  diagram  the  thermal  capacity  of  the  semiconductor  is 
disregarded,  since  it  is  much  lower  than  the  thermal  capacity  of  the  heat 
sink  and  is  therefore  significant  only  in  the  case  of  very  small  pulse 
times . 

The  thermal  resistance  from  junction  to  heat  sink  includes  both  the 
thermal  resistance  from  junction  to  case  and  the  thermal  contact  resistance 
to  the  heat  sink.  Provided  a heat-conducting  paste  is  used  between  semi- 
conductor and  heat  sink,  experience  shows  that  Rj_hs  may  be  set  equal  to 
l*hs— amb»  since  in  the  overwhelming  majority  of  cases  there  will  be  the  same 
thermal  resistance  for  the  heat  sink  selected  as  for  the  semiconductor  in 
question. 

When  a load  is  applied,  there  will  be  a very  rapid  heating  up  of  the 
semiconductor,  which  will  give  rise  to  heating  up  of  the  heat  sink  with  a 
time  constant  corresponding  to  the  product  of  the  thermal  resistance  of  the 
heat  sink  and  its  thermal  capacity. 

The  temperature  will  approach  its  maximum  value  more  or  less  closely 
depending  on  the  relation  between  time  constant  and  pulse  width. 
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When  the  power  is  cut  off,  the 
junction  will  very  rapidly  approach 
the  temperature  of  the  heat  sink, 
and  after  that  its  temperature  will 
decline  with  the  aforementioned 
time  constant  for  the  heat  sink. 

How  far  the  temperature  will  drop 
depends  on  the  ratio  between  the 
"off"  time  and  the  thermal  time  con- 
stant . 

Figure  20  shows  a theoretical 
temperature  variation  for  the  junc- 
tion . 

From  what  has  been  said  here, 
the  load  applied  can  be  calculated 
as  a function  of  the  pulse  time  T, 

and  with  the  junction  temperature,  the  two  thermal  resistances,  and  the 
heat  capacity  of  the  heat  sink  as  parameters.  If  the  junction  temperature 
is  set  at  the  maximum  permitted  (in  the  calculations  that  follow,  150°C) , 
and  with  a defined  maximum  ambient  temperature  (40°C) , the  maximum  permis- 
sible load  can  be  calculated,  since  Rhs-amb  and  chs  are  determined  by  the 
heat  sink  chosen  as  a function  of  the  pulse  time  T and  with  the  "on-off" 
duty  cycle  as  primary  parameter,  according  to  the  formula: 


Figure  20.  Temperature  Variation 
at  the  Junction. 


AT. 


i max 


Rj-hs+ 


1 -e 

-T/$hja-aJnb^  * 


chs 


A derivation  of  the  expression  is  given  in  Appendix  2. 

Figure  21  shows  the  relation  between  the  maximum  permissible  load  in 
constant  service  (P^q)  and  the  maximum  permissible  load  with  a specified 
pulse  time  (P)  as  a function  of  the  pulse  time  (T)  in  relation  to  the  given 
time  constant  for  the  heat  sink  (x)  and  with  the  pulse  duty  cycle  t\/7  as 
a parameter  --  P/Pdc  = f(T/x,  tj/T) . The  time  constant  t of  the  heat  sink 
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is  equal  to  the  product  Rhs-amb  * Chs . The  curves  are  calculated  and  are 
valid  only  with  the  stipulation  that  the  thermal  resistance  from  heat  sink 
to  ambience  is  equal  to  that  from  junction  to  heat  sink. 

All  curves  give  a temperature  rise  from  junction  to  ambience  of  110°C. 

Since  the  junction  temperature  is  difficult  to  measure,  in  order  to 
verify  the  theory  the  same  curves  are  plotted  in  Figure  22  with  the  thermal 
resistance  from  the  junction  to  the  heat  sink  left  out.  Naturally  these 
curves  give  considerably  higher  permissible  loads,  since  the  momentary  rise 
in  junction  temperature  in  relation  to  the  temperature  of  the  heat  sink 
when  the  load  is  applied  is  left  out.  The  curves  show  a temperature  rise 
from  the  heat  sink  (the  hottest  point)  to  the  ambience  of  55°C. 

Measurements  were  done  on  this  basis  for  a number  of  heat  sinks,  and 
various  points  on  the  curves  were  rechecked.  It  was  found  that  all  measure- 
ments gave  temperatures  within  ±2°C  of  the  permissible  temperatures.  Fig- 
ures 21  and  22  are  shown  on  the  next  two  pages. 


Figure  21.  The  curves  show  the  maximum  permissible  load  in  intermittent  service  in  relation  to  the 
maximum  permissible  load  in  constant  service  as  a function  of  the  ratio  of  pulse  time  to  the  time  con- 
stant of  the  heat  sink.  All  curves  give  a temperature  rise  of  110°C  at  the  junction.  [The  subscripts 
in  the  legend  signify:  J - junction,  CMC-  - ambience,  K - heat  sink,  DC  - continuous  service,] 
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10.  TEN  EFFECTIVE  WAYS  OF  DECREASING  THE  LIFETIME  OF  A POWER  COMPONENT 

—^i  » ■ .i  ■n.i  — —ii——  i j — — »w wip—  o—  — — ■ ^ i »wi wmmem 

1 . Forget  that  an  increase  in  load  produces  the  greatest  temperature 
rise  across  the  greatest  thermal  resistance  (possibly  not  the  heat  sink) . 

2.  Forget  the  thermal  contact  resistance  between  semiconductor  and 
heat  sink. 

3.  Dispense  with  the  use  of  a thermal  compound, 

4.  Cover  the  heat  sink  on  the  top  side. 

F,  Forget  the  possibility  that  the  apparatus  will  be  '’packed  down" 
while  in  service. 

6.  Mount  the  heat  sink  horizontally  or  in  such  a way  that  the  ribs  are 
parallel  to  the  horizontal. 

7.  Use  a bright -surfaced  heat  sink  and  calculate  with  the  data  sheet's 
curves  for  a black  anodized  one. 

8.  Forget  the  possibility  of  solar  radiation. 

9.  Forget  the  radiation  from  other  hot  components  or  surfaces. 

10.  Calculate  with  a maximum  ambient  temperature  of  25°C,  especially 
when  the  heat  sink  is  mounted  in  a cabinet. 


APPENDIX  1,  COMPUTER  CALCULATION 
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The  input  program  discussed  in  Section  8,4  is  shown  on  this  and  the 
next  two  pages.  The  program  is  in  "Basic."  Lines  90  to  620  and  1300  and 
1310  are  input  lines.  Lines  620  to  960  are  calculation  of  the  necessary 
constants  and  the  rest  are  output  lines.  An  example  of  such  an  output  is 
shown  on  page  SO.  The  first  symbol  after  the  line  number  indicates  whether 
an  impedance,  admittance,  capacity,  or  current  is  involved;  the  next  two 
figures  indicate  the  nodes  between  which  the  admittance,  e.g.,  is  situated; 
the  last  number  shows  the  size  of  the  component.  This  output  can  be  used 
directly  as  input  to  a circuitry  program. 

COL  1 08:40  09/30/74 

90  FILES  DCFIL 

100  PRINT  " THIS  IS  A PROGR.AM  FOR  CALCULATION  OF" 

110  PRINT"  TEMP.  DISTRIBUTION  ON  A HEAT  SINK 

120  PRINT 

130  PRINT  "WHAT  IS  THE  PROFILE  LENGTH  - L,  RIB  HEIGHT  - H," 

140  PRINT  "RIB  INTERVAL  - A,  DISTANCE  FROM  CENTER  LINE" 

150  PRINT  "TO  FIRST  RIB  - C,  THICKNESS  OF  RIB  1ST  AT" 

160  PRINT  "THE  BASE  - SI , AND  2ND  AT  THE  OUTER  END  - S2,  MEAS . IN  METERS." 
170  INPUT  H3,H4,H1 ,H2 ,S1 ,S2 

180  PRINT  "HOW  MANY  RIBS  ARE  THERE  ON  EACH  HALF  OF  THE  HEAT  SINK"; 

190  PRINT 
200  INPUT  N 

210  PRINT  "WHAT  IS  HALF  THE  THICKNESS  OF  THE  BODY" 

230  PRINT  "BETWEEN  CENTER  AND  1ST  RIB,  BETWEEN  1ST  AND  2ND  ETC" 

235  PRINT  "VALUES  ARE  PUNCHED  IN  ONE  AT  A TIME" 

240  FOR  K=1  TO  N 
250  INPUT  M(K) 

270  NEXT  K 

280  E=INT (H3/0.03+0. 5) 

290  F=INT (H4/0.01+0.5) 

300  PRINT  "LENGTH  IS  DIVIDED  INTO  I";E;"NODES" 

310  PRINT  "RIBS  ARE  DIVIDED  INTO  I";F;"N0DES" 

320  PRINT  "IS  THIS  OK  YES  OR  NO"; 

330  INPUT  XS 

340  IF  XS-"YES"  THEN  380 

350  PRINT  "WHAT  SHALL  THE  NUMBER  OF  NODES  BE  IN" 

360  PRINT  "RELATION  TO  LENGTH  AND  RIBS"; 

370  INPUT  E,F 

380  PRINT  "WHAT  IS  THE  AMBIENT  TEMPERATURE," 

390  PRINT  "WHAT  IS  THE  TOP  TEMPERATURE  OF  THE  HEAT  SINK," 

400  PRINT  "AND  WHAT  IS  THE  TEMPERATURE  GRADIENT  BETWEEN  THE  RIBS." 

410  INPUT  T 3 , T 1 , T 2 
415  T3=T3+273 

420  01=(T3+T1)“3+ (T3+T1) *T3~2+T3* (T3+T1 )~2+T3  3 

430  02=  (Tl+T2+T3}~3+  (T1+T2+T3)*  (T3+T1)*‘2+  (T3+T1)  * (T1+T2+T3  2+ O'3+Tl  T 3 


A? 


443  pniriT  "IS  THE  HEAT  SINK  MOUNTED  VERTICALLY  OR  HORIZONTALLY"; 

459  INPUT  XS 

«to  "VERTICALLY"  THEN  510 

470  X 5«2 
4*3  VI  *9. 71 
4)0  L7*H4/(F*“ I 3 
5,10  GO  TO  549 
St 0 U-H3/S 
579  )'t  »3  *55 
530  L2*L1 

549  V2-9.55 

550  P»1MT  "WHAT  IS  THE  RADIATION  NUMBER"; 

560  INPUT  tl 

570  T2«l I .34*tl /<5. 77-2*51 > 

550  Z|«Et*lE-8 
540  22»J2»i£-S 

MO  ppint  "WHAT  IS  THE  THERMAL  CONDUCTIVITY  OF  THE  VATERIAL"; 

610  l .•‘PUT  G 

670  Xl«2*G*H3/E 

630  X2«G*<E-I >/H3 

640  X3»G»H4»CE-1  >/<F-n/HJ 

650  X«-G»H3»<F-l J/H4/E 

669  X6*2*5*H3#,’I  /E 

670  X7»3*V2*H4#H3/E/<F-t > 

660  ra«£j •H3»Ot /E 

719  r<4)«Zt*H4#K3*01/tt»{r-l>» 

720  Y<8>»Z2*K4*:i3,‘03/CE*tF-t>> 

733  Y<4)«2#G#Ht 1 >*H3/H2/E 

740  YC  t-0  )»G*H2*M<  t >*<E-t  >/(2*H3> 

760  A«F*E 

770  ?■ A*E 

780  SCPATCH4I 

790  FOP  )>t  TO  It 

830  0»E*<F*l > • < K - t »♦£ 

820  IT  I>1  THEN!  856 

830  H5-H2 

540  GO  TO  360 

350  H5-KI 

360  FOP  I«l  TO  E 

870  FOP  J*l  TO  <F*I > 

e?o  if  j>2  then  939 
890  xc i *2>«xt • ut::>/H5 
695  IF  :>1  THEN  90  5 
930  XC  I » I )«X2»:i(X>#lt5/2 
902  30  TO  910 
905  XI  t * I )«X2»:i(”)*H5 

919  x<?*  t >«x2*.i(:t>#st 

920  GC  TO  940 

930  r(J.|)»X3,(S2*<Sl-S2>«<l-tJ-?>/F>> 

940  F J<7  the:*  960 

959  ~<J»2)*X4*<  S3*<  St -S2>«<  t-<J-';  >/.  J> 

960  IF  l«t  THEN  980 

970  '’PI  TE*  I , *Y*  J '*•♦  ( t-t  > ♦<  J - 1 ) *EJ  2*<  I - I > ►(J'l  > •£♦ IIX( J*  I >1 
9Cd  IF  .>•  < F*  t ) THIN  I 140 

990  ’'•’iTi#i.*Y-;r*i*tj-i)«En*i*J*Eix<j,2)» 

1000  IF  J« 1 THEN  1070 

iota  D«7»l *E*< J-l  ) *E 

1020  Y(3>«X7«S'3P(50?.CTI/L2>) 

1039  V’P.l  TE*  I , *Y";  D!  0 I Y(  3 > 1 

10  40  IF  H*N  T’-i  IN  1 C 60 

1050  T'UTEf  l,*Y-JDJ0JYC4>J 

1060  GO  T0  1139 

1073  IS  XS«»2  THEN  1090 

1030  L)  •IH5,H3/E>/<H3/E*1I5» 

i 0 9 0 Y(5)*X6*.i5*S3X(S3~(TI/l.l  ) > 

1193  Y<6>»X0*H5 

1110  "PITE4 I . -Y* J n*I ; 0JY( 5>) 
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1120  WRITE#1  , •,Y*';a*IJ0m6>J 

1130  NEXT  J 

> 1 40  NEXT  I 

1150  FOR  J-!  TO  E 

1160  IF  K«1  THEN  ? 19S 

1170  B»2*E+J*P*<K-2> 

1180  YC7J-X1 •K<K)/Ht 

1190  write# i,“Y-jb;b*a;yc7>j 

1195  IF  K-N  THEN  1280 
1200  FOR  1-1  TO  <F-1> 

1210  C-3*E*<K-1 >*P*JfE*CI -1 > 

1220  vrite#i,-y-;ciop;y(8>; 

1230  NEXT  1 

1240  IF  K>l  THEN  1280 

1250  VRITE#1 > “Y-i Ji J-EiY(9>I 

1260  IF  J-E  THEN  1290 

1270  WRITE# 1,“Y“I JiJ+i;YC 10)1 

1280  NEXT  J 

1290  NEXT  K 

1300  PRINT  "HOW  GREAT  IS  THE  POWER  IN  WAITS  (1/8V 

1310  INPUT  I 

1320  urite#i,-i-;o;ui/2; 

1325  VRITE#!  4 “I -;o  J2;  I/W 
1330  RESTORE# 1 
1340  END 
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APPENDIX  2.  DERIVATION  OF  THE  EXPRESSION  P * 


AT . 



-ti/Rhs-amb  * chs 


Rj  -he  + Rhs -amb'l-g_T/Hh3-eun'b  • ^ he 


In  the  following  derivation  of  the  above  expression  the  thermal/elee- 
trical  equivalents  shown  in  the  table  below  are  used. 


Thermal  Magnitude 
P power 
T temperature 
R thermal  resistance 
C thermal  capacity 


Electrical  Analogy 

v I current  A 

’C  'v  V voltage  V 

JC/W  <v  R electrical  resistance  H 

J/°C  *v  C electrical  capacity  F 


'this  means  that  in  the  present  case  (cf.  Figure  IS)  we  have  the  follow- 
ing analogy: 


R1  R J-hs 

^2  R hr-amb 

C % C ns 

The  figure  below  gives  the  analogous  electrical  equivalent  diagram. 


r*  Hi 


r.  v‘j 


Figure  23.  Analogous  Electrical  Equivalent  Diagram  for  Intermittent  (Pulse) 
Service  of  the  Heat  Sink. 

Since  we  are  only  interested  in  the  maximum  temperature  rise  in  the 
junction,  in  the  derivation  of  the  expression  we  shall  concentrate  only  on 
V'max  and  Vmax  respectively  and  not  v'  and  v. 


I 


• f wwwwww ^T? w* ’"'V v* ■'  ,"n  vr,? gniy .^» Hfp!|||^ 


Figure  24.  Analogous  Maximum  Load  for  Heat  Sink 
as  a Function  of  the  Time.  [For  OP  read  TOP; 
for  N."D  read  BOTTOM.] 


Since  what  we  hkve 
here  is  a current 
generator,  the  voltage 
at  R;  will  contribute 
with  a constant  voltage 
drop  in  the  period  tj 
and  no  voltage  drop  in 
the  period  t2 . 

To  facilitate  the 
calculations,  Vmax  is 
first  found  instead  of 
V 'max ; i.e.,  we  get 
Vmax  expressed  in  terns 
of  i,  t ] , and  t2 . Then 
the  voltage  over  R; 
can  be  added  to  the  re- 
sult to  get  V'max. 

Lastly,  from  this 


expression  i can  be  found  as  a function  of  V'maX>  Rj  and  R2 , C and  t3  and  t2 , 
ind  these  can  be  converted  back  te  the  thermal  magnitudes. 


Figure  24  shows  the  voltage  v as  a function  of  the  time, 
The  calculation  is  shown  below, 

The  first  voltage  maximum  can  be  expressed  as 

v», 


'max 


**  i*Ri  (1-e  , 


The  first  minimum 


V‘mln  = V‘  -*e 


. -t,/R,C 


‘max 


The  second  maximum 

V. - * Vlmln 


‘max 
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wm 


The  second  minimum 


v,  * v2  .S~t*/R»c 

i«<n  ’max 


The  third  maximum 


V*max 

The  third  minimum 

The  fourth  maximum 

V»  „ 
max 


V,  . « V|  .e“t,/R,C 

min  ’max  » 


i*Rj (l-e"t,/R,C)  + Vlmine”t,/R*C  * 


etc . 


Substituting  V3  for  V4  , V2  for  V3 , etc.,  we  get 

Vkmax  “ i*RaU-*“t,/R,C>-Mi*Ri  (l-e“t,/R*C)+(i*R»  (l-e“t,/R*c)  + 

i.RIU-e't‘/RlC).e‘tl/RlC.e-ta/R*c)*“t>/R*c.*‘t*/R*c)- 

e-tj/RjC.e-t2/RjC 

Reducing  this  expression,  we  get 

i*R*  (l-e"t*^R,C)  t 4 


V. 


max 


+ e-3(ti+t2)/RjCj 


Since  the  expression  in  sauare  brackets  is  a quotient  series,  Vmnv 
be  found  by  allowing  the  series  to  approach  the  infinite. 

l-e‘tl/RjC 

Vmax  i#  R ”7"e” T^T+t*7 ' 

The  voltage  field  across  Rj  is  added. 

,_„-tj/R*C 

v.  - 1 ,r,+r2  ■i^-rr~TTr’TT— f*  • 

j max  j it j+t2 ) KjC 


can 


Since  we  wish  to  find  i as  a function  of  tj  and  t2 , we  divide  the 
parenthesis  through  and  we  have 


jmax  ______ 

R,4R* 


The  electrical  magnitudes  are  replaced  with  the  equivalent  thermal 
magnitudes,  and  we  have  the  expression 


^Tjmax 

r + » l-e~tiyRha»am^Sag 

J— he  rie  — ajuT^^-T/Rhe-am)?  C hs 
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